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^—^ Abstract 

C^l ■ We investigate correlation functions in a periodic box-ball system. 

For the two point functions of short distance, we give explicit formulae 
obtained by combinatorial methods. We give expressions for general N- 
' point functions in terms of ultradiscrete theta functions. 

1 Introduction 



> 



Quantum integrable systems such as quantum integrable spin chains and solv- 
able lattice models are systems whose Hamiltonians or transfer matrices can 
be diagonalised and for which eigenstates or free energies can be explicitly ob- 
' tained pQ. To investigate physical properties of these systems, such as e.g. the 

^\ . linear response to external forces, however, we further need to evaluate correla- 

tion functions for these systems. This is one of the main problems in the field of 
quantum integrable systems and in fact, obtaining correlation functions is even 
fairly difficult for the celebrated XXZ model or the 6 vertex model [3J. 
ON . A periodic box-ball system (PBBS) is a soliton cellular automaton obtained 

by ultradiscretizing the KdV equation [31 0] . It can also be obtained at the q — > 
limit of the generalized 6 vertex model [3 [6] . Hence, from the view point of 
quantum integrable lattice models, it is interesting and may actually give some 
new insights into the correlation functions of the vertex models themselves, to 
obtain correlation functions of the PBBS. In this paper, we give expressions for 
A-point functions for the PBBS, using combinatorial methods and the solution 
for the PBBS expressed in terms of the ultradiscrete theta functions. 

The PBBS can be defined in the following way. Let L > 3 and let Hl = 
{/!/:[£]- {0,1} such that ^({l}) < L/2 } where [L] = {1, 2, . . . , Z}. 
When / 6 SIl is represented as a sequence of 0s and Is, we write 

/(l)/(2).../(L). 

The mapping Tl '■ — *• is defined as follows (see Fig.[lJ: 



o 



X 
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f= 0011001111000111000 



T L f= 1000110000111000111 



Figure 1: Definition of Tl for / e FIl 



1. In the sequence / find a pair of positions n and n + 1 such that f(n) = 1 
and f(n + 1) = 0, and mark them; repeat the same procedure until all 
such pairs are marked. Note that we always use the convention that the 
position n is defined in [L], i.e. n + L = n. 

2. Skipping the marked positions we get a subsequence of /; for this sub- 
sequence repeat the same process of marking positions, so that we get 
another marked subsequence. 

3. Repeat part 2 until one obtains a subsequence consisting only of 0s. A 
typical situation is depicted in Fig. [T] After these preparatory processes, 
change all values at the marked positions simultaneously; One thus obtains 
the sequence T^f. 

Sometimes we shall write T l L f for T L (- ■ ■ {T L {T L /))). The pair (fl L ,T L ) is called 

t 

a PBBS of length L [H[7]- An element of £Il is called a state, and the mapping 
Tl the time evolution. 

An iV-point function of the PBBS with M balls may be defined as follows. 

( Sl ,s 2 ,...,s N ):=^- Y, ^ H ^f( Sl )f( S2 )---f(s N ) 

where Q L;M := {fefl L \ (^({l}) = M }, Z H := E/ e n £iM eZ^^f) and 
Hk(f) is the kth energy of the state /, which is proportional to the number of 
fcth arc lines defined when determining the time evolution rule [4], or the fcth 
conserved quantity of the PBBS [S]. (Note that Hk(f) is essentially equal to the 
energy function for the transfer matrix of the crystal lattice models with fc + 1 
states on a vertical link [3 [7].) Noticing the fact that Hl-,m = \_\y ^y, 

( S1 ,S 2 ,...,S N ) = ^-Y1 E ^ LlHkU1 f(si)f(s2)---f( S N), 

where Y are partitions of M corresponding to the conserved quantities of the 
PBBS. (See Section H) Since, for /< e fl Yi (i = 1,2), v fc H k (fx) = H k (f 2 ) 
(k = 1,2,3,.. .) implies Y\ — Y 2 and vice versa, by choosing a state fy in fly 
we can write 
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Thus, to obtain correlation functions of PBBS, we have only to evaluate those 
on the set fly' 

(s 1 ,s 2 ,... ! s^)r := —r V f(si)f(s 2 )---f(s N ). (1) 
|1V| f&lr 

We also point out that if we put v fc, v /, Hf-(f) = 0, TV-point functions become 
trivial; 

. l-nCm-n M(M-l)---(M-N + l) 

( Sl ,s 2 ,...,s N )= = L(L _ 1) ... {L _ N + 1) ■ 

In the following sections we shall evaluate fl}. 

First we summarize some useful properties of the PBBS. We say that / has 
(or that there is) a 10-wall at position n if f(n — 1) = 1 and f(n) = 0. Let the 
number of the 10- walls be s and the positions be denoted by a\ > ai > ■ ■ ■ > a s . 
Then, we have the following proposition: 

Proposition 1 (|9j) 



(Tlf)(n) = 



'ln+1 'ln+1 'In ' '/n> 



max 

<6M 



i-1 



Zi 



+ ^2min{W i ,W 7 } + W i 

3=1 

min{W 4 ,V^}, 



J2 miE 

1=1 3=1 

2 ' 



(2) 
(3) 



Z< = i-^2min{W l ,M^}, 

3=1 

where Wi denotes the amplitude of the "soliton" corresponding to obtained 
by the procedure explained in 0. 

The set {Wi}*_ 1 consists of quantities of the PBBS and rf n is the ultradis- 
crete theta function [TU] . We shall use Proposition [T] for determining iV-point 
functions in Section [3] 

Next we introduce two procedures which are important in this paper. For a 
given / £ fi^, a state Ef = E(f) is defined to be 

f(n + 2k) 



(Ef)H = 



(1 < n < a s - 2), 

«s-fe+i — 2fc + 1 < ?i < a s _fc — 2k — 2 
(fc = l,2,...,s-l) 
f(n + 2s) (ai —2s<n<L — 2s), 



(1 < n < a s _i - 3), 

a s _fc — 2k < n < a s _/j_i — 2k — 3 



/(n + 1) 

^ + 2fc+1 ) V (fc=l,2,..., S -2) 
/(n + 2s - 1) (ai - 2s + 1 < n < L - 2s) 



(a s > 1) 



(a. = 1) 
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The mapping E : £1^ — > 0^_2s is called the 10-climination. Ef is a subsequence 
of / obtained by eliminating all 10-walls in / simultaneously. For example, 

/ = 0011111000011111000000001111001100001111100111000011110000000000, 
Ef = 001111 0001111 0000000111 01 0001111 Oil 000111 000000000 
= 00111100011110000000111010001111011000111000000000. 

Its inverse process is called the 10-insertion, 1(31,32, • ■ • ,jd) = h°h{ji,j2, ■ ■ ■ ,3d) ■ 
CIl — ► &L+2(d+s) where s is the number of 10-walls in / e SIl- The 10- 
insertion is defined as follows: Shifting the origin if necessary, we can assume 
that f(L) = 0. For {ji, j 2 , ■ ■ ■ , jd} (1 < h < h < • ■ • < 3d < L + d), the 

^L+2d is defined as 



mapping h(ji,h, ■■■,3d) ■ 

(h(ji,32, ■ ■ -,jd)f)(n) 

1 (n = L + 2d- j k - k + 1), 

(n = L + 2d-j k -k + 2), 

/(») {l<n<L + d-j d ), 

f(n - 2(d - k + 1)) (L + 2d - j k - k + 3 < n < L + 2d - j k -i - k + 1), 

f(n -2d) (L + 2d - jx + 2 < n < L + 2d) 

where fee [d]; furthermore, Ii : &L+2d ~^ ^L+2(d+s) is defined to be 

' 1 (n = g k + 2(s-k) + l), 

(n = g k + 2(s-k) + 2), 

{h.f')(n) = { /'(«) (l<n< fla ), 

f'(n — 2(s — k + 1)) (g k + 2(s - k) + 3 < n < g k -i + 2(s — k) + 2), 

f'(n-2s) (g 1 + 2s-2<n<L + 2{d + s)) 

where k 6 [s], f = h{ji,j 2 , . . -,jd)f G ^L+2d and 

= max | m € [L + d] m = at-l + f{r£ [d] \ L + d — j r + 1 < m } j, 
Sfc - <7fc + tt{ r G [d] | L + d - j r + 1 < g' k }• (4) 
For example, 

/ = 0011100111000001101000111000000, 
/i(3,ll,25)/ = 001110011* 1000001101000* 1110000* 00 
= 0011100111010000011010001011100001000, 
7(3,11,25)/ = 001110001110100000011001000010111000001000 



where 10 and denote the inserted 10 at / ^ h{j\,j 2 , ■ ■ ■ ,jd)f and h(ji,j 2 , ■ ■ ■ ,jd)f 
h{h(.h,h, ■ ■ ■ ,3d)f) respectively. 
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2 One and two point functions obtained by com- 
binatorial methods 

We assume that Y denoting the conserved quantities of / £ f2y , is the partition 
(Pi,Pw- ,Pi,P 2 ,P 2 ,--- ,Fi,...Pt,Pt,-~ ,Pt) 

V v ' S v ' S v ' 

n\ ri2 ri£ 

where P% > P 2 > ■ ■ ■ > Pg > 1. Note that Y is a partition of M, i.e. M — 
S»=i n iPi- As mentioned in Section [fl we consider TV-point functions Q of the 
PBBS, 

(s 1 ,s 2 , ■ ■ ■ ,s N ) Y = |^— | 51 f(si)f(s 2 )---f(s N ). 

The value of fiy is already known: 
Proposition 2 ([II]) 

\ = J± ( L + m-l\f L x + n 2 -\\ { Li-x + m-1 
' Y ' L V n i / \ "2 / \ 
where L„ = L - 2M, L { = L Q + Y.)=i 2n j( p j ~ p i+i) and P e+1 = 0. 

Since the TV-point function (si, Si + di, . . . , Si + dAr-i}y does not depend on 
the specific site S\ (because of translational symmetry), we denote 

C Y (di,d 2 , . . . ,djv-i) = (si, si + di, . . . ,si + djv-i)y 

where 1 < d\ < d -2 < ■ ■ ■ < djv-i < Note that Cy(0) denotes the 1-point 
function (si)y. 

Proposition 3 

M 

Cy(0) = -. 

Proof Since £^ =1 f( n ) = M > 

r 1 /en y ti=i 



□ 



Next we consider the 2-point functions. 
Proposition 4 

M-s 



Cy(l) 

where s = £) i=1 
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Proof Since £^ =1 f(n)f(n + 1) = M - s, 



1 L 

LC Y (1) = — E /( n )/( n + !) = M - 3 - 
' Y ' fen Y n=i 



In order to investigate Cy(2), let us put 
&i := 



We also define 



I otherwise, 
Pi 

L:=L-2k l (=L-2s), 
Pi 

N Y (2) :=^ki{i-2). 

i=3 



V fa := { / e n Y I £/ = /o } , 
G 2 (/) :=tt{nG[L]|/(n)/(n + 2) = l}. 



The following lemma is the key to evaluating Cy(2). 

Lemma 1 

Let 

V%> :={feV fo \G 2 (f) = N Y (2)+j}. 
Then, ifV fo + 4>, V fo = Ujio V h and 



y(fci-i) 
Jo 



fci! 



where 



\ /fei+J-l N 
\i=0 /V i=0 y 



X i 



e n 



i 



o<n<...< s .< fcl +j-i h= i {2k 2 +i h )(2k 2 +i h + 1) / 

*h + K*h+l 



□ 



(5) 



Proof When / G V/ , there exists a set of positive numbers {ji}^ (1 < ji < 
32 < ... < jki <L + k\) such that 

/ = HhJ2, ■ ■ ■ ,ifei)/o- 
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By examining the positions of 101 and 111, we find that 

G 2 (f) =AV(2)+ 7 + J 

where 7 = j(fo; {jijfei) i s the number of 10s inserted into the positions adja- 
cent to consecutive Is, and J = |j{ i € [d — 1] \ji + 1 = ji+i }. (See the table 
below.) For example, 

f = 001110000100110000 
and / = 7(5, 6, 14, 15, 18)/ , then 

/ = 0011110100010100011000111010100000 
( = OOlll|l^OOlOlOOO10OOll|l^lOlOOOOO ). 

In this example, k\ = 5, k 2 = 3, Ny(2) — 3, 7 = 2 and J — 2. Since 
< 7 + J < ki , we have the decomposition Vf = U jlo ^0 ' 



fo 


00111000 


G2U0) 


1 


f = I(fc)/o 


00101110000 

(fc = 7) 


00110110000 

(fc = 6) 


OOlll|01OOOO 

(fc = 4) 


00111000100 

(k = 2) 


Ga(/) 


3 

(7 = 1, J = 0) 


2 

(7 = 0, J = 0) 


3 

( 7 =1, J = 0) 


2 

(7 - 0, J - 0) 



To know V^ o , we have only to count the number of states with 7 + J = j. 

For hi = 1, V/J = i. Since there are ki sets of consecutive Is, Iki states 
have 7 + J= l (7 = 1, J = 0) and the other i — 2k 2 states have 7 + J 

(7 = 0, J = 0). 

For fci = 2, let / = I{ji 1 j2)fo- As was seen in case k\ = 1, there are 2fc 2 
positions at which 7 + J can be increased by one. If one 10 pair is inserted in one 
of these positions, then there are 2k,2 + 1 positions for the other pair to increase 
7+ J by one, and L—2ki positions not to increase it. On the other hand, if one 10 
pair is inserted at one of the L—2ki non-increasing positions, then there are 2^2 + 
2 positions for the other pair to increase 7+ J by one, and L— 2k 2 — 1 positions not 
to increase it. Hence, considering duplication of insertion, there are (2fc 2 )(2fc 2 + 

l)/2! states with 7 + J = 2, (2k 2 )(L - 2k 2 ) + (L - 2k 2 )(2k 2 + 2) /2! states 

with 7 + J = 1, and (L - 2k 2 )(L - 2k 2 - l)/2! states with 7 + J = 0. 

In general, we can proceed in a similar manner and, referring to the chart 
in Fig. we obtain (|5]). □ 







Proposition 5 



C Y {2) 



ki / Pi 

I>j * (!>(*- 2) + (k-j) 

j=0 \i=3 

ki 

3=0 



7 



2k, +4 2 
2t, +5 2k, +3 2 

2i, +6 2k, + 4 2/t, 



2*, +2 
L-2k, 



Figure 2: A chart corresponding to 7 + J in the proof of Lemma [TJ 



Proof From Lemma [TJ we see that if Vf ^ 



E E /(")/(" + 2 ) = E S ( iVy(2) + (fcl ~ j) ) 



/eV) n=l 



and 



j'=o 



Since the right hand side of the last equation does not depend on fg, and since 
any state / G fly belongs to some Vf , we obtain 



£CV(2) = J-r E E/W/(n + 2) 



E^(^y(2) + (fc 



1 - J. 



/Efiy n=l 



fel 

E^ 

j=0 



□ 



For Cy{d) (d > 3) we can use similar arguments based on elementary com- 
binatorics. However, the expressions become more and more complicated when 
the difference d increases. Instead in the next section we shall use Proposition!]] 
to obtain expressions for general iV-point functions. 
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3 TV-point correlation functions for the PBBS 

Let the state f and the set X Y C Z™ 1 x Z'j 2 x • • • x Z 1 ^ (= 7L\) be 



and 



fo = 000^_00, 

L 



1 < Xi(l) < Xi(2) <■■■ < Xi(n,i) < 7 4 _i + n, 

(i = i,2,...,e) 



(6) 

We define the state fx recursively as 

f r .= J(0)---/(0) AX-)U-x (j = 1,2,. ..,£), 

V v ' 

-Pf-j + l--P<!-j + 2-l 

/x := ft 

where Xj = {^-(fc)}?^ C X E Xy . Note that, from the definition of an 
10- insertion, 7(0) is the procedure needed to insert s between 10: 

/ = 0011100111000001101000111000000, 
7(0)/ = 001110001110000001100100001110000000 

and /j £ !2y by construction. We also define f2y by 

Q Y := { f x | X e *y } • 

Lemma 2 



1 L 

( Sl , S2 ,..., s ^) y = -_- ^ 2/(fc + «i)/(fc + s 2 )"-/(fc + «iv). (7) 

Proof By virtue of the definition of fx , &y is the set of states with conserved 
quantities Y and the last entry of the 10 sequence is one of the 0s that are not 
marked in the time evolution rule, i.e., fx(L) = (1l/x)(7) = 0. By defining the 
shift operator S by (5/)(n) := f(n + 1), and {S k f) := S(S k -\f) (k = 1,2,...) 
with 5°/ := / and for sets 

S k n Y := { S k f x | X 6 X Y } (k = 1, 2, . . . , L), 

we find 

v f eO y , tf{fc|/eS fe Qy (k = l,2,...,L)}=L . 
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Note that S L f = f. Since \Q Y \ = t^Py\, 

(s 1 ,s 2 ,...,s n ) y = |^-| /( s i)/( s 2)--'/( 



Thus we obtain ©. 



/en* 



Lin 



jnE E M)/(*a)---/(«jv) 

1 fc=1 fes k n Y 
1 L 

^EE +*)/(** + *)•••/(** + *)• 



| * 1 fc=i / e a 



□ 



Proposition 6 

For X £ Afy > fx is explicitly given as 



where 



f x (n)=u° n (X) 



<(X) := V t - + \(X)- V t n+1 (X)-r 1 t - 1 (X)+r 1 t n (X), 



max 

ie[«]i je[ni] L i=l fc=l 



Y E m «(* p i - n - x ^ k ) + L + k + 1 + -y) 



^ E E E E m ik^ik]h-m jh 
i=l k=l j=l h=l 



(8) 



Zi 



2 fiij^kh ~\~ Pmax[i,j]i 

(i i 
^ E" ' • E 
i=l j=i+l 



Proof From Proposition [TJ /x is determined by the parameters W n and a n 
(n = 1,2, ...,s). Here W n is the amplitude of the nth soliton and a n is its 
position, i.e. the position of the nth 10- wall, counting from the right. From the 
definition of the position and of the amplitude of a soliton, it follows that both 
can be determined from 10 insertions. Because of the way fx was constructed, 
the set {xj(k)} f J_ 1 corresponds to the position of rij solitons with amplitude 
Pj, though it does not directly gives their position. Hereafter we shall refer to 
a soliton with amplitude P as a P-soliton. By considering the relation between 
the position of a soliton and 10-insertions, we find that the position of the kth 
Pj-soliton counting from the right is L — xf' (k) + 2, where x"p (k) is determined 
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recursively: we define Xj(k) (i £ [£], j £ [i], k £ [rij]) as 



xf\k) := Xj (k) + (P, - P i+ i)(2/3,(fc) +2k-l) 



where 



2(P s -P l+1 )af(k)-k + l 



s=j+l 



af{k) := r £ [m] | L;_i + m - Xi (r) + 1 > gf (k) }, 



+ rii — x 



(k) + l), 



gf\k) := max < m £ [Lj_i + m] 



B f~ 1] ' 



x: >(k) + i 



+tl{ r £ [m] I + - Xj(r) + 1 < m } 



Note that ^(1) < 2^(2) < • ■ ■ < xf'irij). 

Recalling the fact that tt{ r £ [d] \ L + d — j r + 1 < .9^. } in (U]) is the number 
of inserted 10s, on the left of the fcth soliton (here we do not count the inserted 
10 s as solitons), the concrete meaning of these variables becomes clear: a.j (k) 
denotes the number of Pi-solitons on the right of the fcth Pj-soliton, and f3j(k) 
denotes the number of solitons with amplitudes less than Pj , to the right of the 
fcth Pj-soliton. 



(0, 



Since {L — Xj (fc) + 2Yj' = i k=\ ^ s * ne complete set of positions of the solitons, 
there exists a one to one mapping p : { (j, fc) | j £ [£], k £ [nj]} — > [s] such 
that 

a P (j,k) = L - xy(k) + 2. 
From these recursion relations we have 

l 

xf{k) = Xj (k) + P j (2/3 j {k) + 2k - 1) + 2p i a fi k ) + l 

i=j+i 

= Xj(k) +2iP j (Pj{k) + (fc - 1)) + p i a f( k ) \+Pj-k+l. 

[ i=3+l J 

Since the position of the fcth Pj-soliton is a, p fj t k)i Wp(i,fc) = p j an( i the set 
of amplitudes of the solitons on the right of the fcth Pj-soliton is nothing but 

{Wh}^' 1 . From the definition of af(k), /3j(k), 



i I, n 



af{k) = ^{w£{W h } P ^ ) - 1 

p j {k) = i{w£{w h y^- 1 



W = Pi), 
W>P 3 ) 



and 

Thus we obtain 



${w£{w h y^ ) - 1 \w = p j ) = i - 1. 



P(j,k)-1 



(fc) = Xj (k) + J2 2 m ™{Wp(j, k ),W h } + W p(j}k) - k + 1. 



h=l 
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Therefore we find a concrete expression of a> p (j,k) , and (JSJ) is immediately ob- 
tained from © and ©. □ 



From Lemma [5] and Proposition [5J we immediately obtain the following 
theorem: 

Theorem 1 

Let X be the set defined in © we have 

L N-l 

C Y (d 1 ,d 2 ,...,d N - 1 ) = — — ^ ^2 U n(X) Y[ u n+di(X), 

I XeXn=l i=l 

for u n (X) = as given in (JSJ). 



4 Concluding remarks 

In this article, we investigated correlation functions for the PBBS and obtained 
explicit forms for 1-point and 2-point functions at short distances. We also 
give expressions in terms of ultradiscrete theta functions for general TV-point 
functions. Investigating their asymptotic properties and to clarify the relation 
to correlation functions for quantum integrable systems are problems that will 
be addressed in the future. 

Finally we should comments on the time averages of quantities in the PBBS. 
The time average: 

, T S L N-l 

Cf(d u <h, . . . , djv-i) = jTzrj ££(T£/)(n) J] (T*/)(n + dj) 

where Tf is the fundamental cycle of / £ depends not only on the con- 
served quantities of the state but, in general, also on the initial state / it- 
self. For example, the conserved quantities of the states fi — 0100100 and 
/2 = 0101000 are the same, but 0^(3) = j and C/ 2 (3) = 0. Hence, in general, 
Cf(di,d 2 , . . . , djv-i) 7^ Cy(di,d2, ■ ■ ■ , rfjv-i) even f° r / S fly. Note that, for 
the 1-point function C/(0), we can easily show that 

v /efV, C/(0) = C Y (0) = ^. 
Acknowledgement 
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A Example of values for the correlation function 

From Theorem [IJ we obtain the following examples, 
(a) L = 12; P x = 3, m = 1; P 2 = 1, n 2 = 2 : 

c 1 i o in c 

Crm = Tr Cy(l) = ~, Cy(2) = -, Cy(3) = -, <V(1,2) = -; 
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(b) L = 14; P x = 2, m = 2; P 2 = 1, n 2 = 2 



q i t; so 

Cy(0) = -, Cy(l) = -, CV(2) = -, C K (3) = — , Cy(l,2) = 0; 

(c) L = 14; Pi =3, m = 1; P 2 = 1, n 2 = 3 : 

3 1 5 69 5 

C(0) = ^ C(l) = -, C,(2) = -, ^(3) = —, ^(1,2) = -; 

(d) L = 14; Pi = 3, m = 1; P 2 = 2, n 2 = 1; P 2 = 1, n 2 = 1 : 
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